Using spectroscopic autoradiography of alpha particles for the quantitative mapping of 226Ra ultra-traces in geo-materials.

The measurement of 226Ra and the identification of 226Ra-bearing minerals are important for studying the behavior of radium in the environment. Various instruments for measuring 226Ra are currently used: among the radiometric techniques that measure in bulk (no spatialization), there are gamma spectrometers and alpha spectrometers. Other instruments such as SEM-EDS can map the chemical elements thus providing information on the distribution of 226Ra, but are limited for ultra-trace analyses on natural geomaterials. Finally, autoradiography techniques can locate radioactivity, but are limited to the identification of the contribution of 226Ra when the 238U series is complete. This study focuses on spectroscopic autoradiography, a method for measuring both the energy of the alpha particle emissions and their positions on the autoradiograph. A gas detector based on a parallel ionization multiplier technology was used for this purpose. Alpha particle energy is dependent on the emitting radionuclides. In order to track the 226Ra, the energy spectrum of the 238U series was studied with modeling software. It appears possible to apply a thresholding on the energy spectrum to discriminate the 226Ra from the first alpha emitters of the 238U decay chain (i.e. 238U, 234U and 230Th, all below 5 MeV). The developed method was applied to a U-mill tailing sample prepared as a thin section. The sample was heterogeneous in terms of radioactivity and was not at secular equilibrium with 238U, as expected. The 226Ra was identified and localized, and different regions of interest were also analyzed with SEM-EDS elements cartography. This revealed 226Ra-rich barite (BaSO4) phases measured at 3 ppmRa on average and containing no uranium; and uranium in siderite (FeCO3), showing a strong 226Ra deficit compared with secular equilibrium. Spectroscopic autoradiography opens up possibilities for the analysis of heterogeneous geological samples containing natural alpha emitters such as 238U and 226Ra: the 226Ra can be localized and quantified at ultra-trace content, and the method developed can also identify newly (young) uranium phases by measuring 238U/226Ra activity disequilibrium.

[18F]-FDG uptake in brain slices prepared from an aged mouse model of Alzheimer's disease using a dynamic autoradiography technique.

OBJECTIVE: 2-[18F]fluoro-2-deoxy-D-glucose positron emission tomography ([18F]-FDG-PET) is a imaging modality that has been used to measure of glucose metabolism in the brain in Alzheimer's disease (AD). Clinically, decreased glucose uptake has been reported in the brain of AD, although the precise underlying mechanisms have not yet been elucidated. To elucidate the mechanisms of decreased [18F]-FDG uptake in the AD by PET, [18F]-FDG uptake in the brain of aged model mouse of AD was investigated using a dynamic autoradiography technique "bioradiography". A X-ray phase-contrast imaging (X-PCI) and a histopathological evaluation were also investigated to elucidate the mechanisms underlying the relationships between decreased [18F]-FDG uptake and the pathological changes in the brain of AD mouse. METHODS: In this study, AD model mouse (5XFAD, APP+/PS1+) were used. [18F]-FDG-bioradiography was conducted in fresh slices of brain tissue under the condition of resting (slices immersed in 5 mM K+ solution) and metabolically active (in 50 mM K+ solution). Amyloid ß42 (Aß42) deposition in the brain of AD mouse was confirmed by X-PCI. In addition, the positive cells of phosphated tau protein (P-tau) and deposition of Aß42 were also examined by immunohistochemical staining. RESULTS: No significant differences were observed between the two groups in the resting condition. In the activate condition of the brain, [18F]-FDG uptake was significantly decreased in AD mice compared to WT mice. In X-PCI showed Aß deposition in the AD mouse, but not in the WT. The AD mouse also showed increased P-tau, accumulation of Aß42, increase in neuronal apoptosis, and decrease in the number of neurons than that of the WT mouse. CONCLUSION: Neuronal damage, and induction of neuronal apoptosis, decreased [18F]-FDG uptake, increased Aß accumulation and P-tau induced neurofibrillary degeneration are observed in AD mouse. In clinical diagnosis, reduction of [18F]-FDG uptake by PET is one of the means of diagnosing the onset of AD. Our results suggest that decreased uptake of [18F]-FDG in the brains of AD may be associated with neuronal dysfunction and cell death in the brain.

Multiplexing Autoradiography.

Autoradiography, the direct imaging of radioactive distribution in tissue sections, is a powerful technique that has several key advantages for the validation of PET radiotracers. Using autoradiography, we can localize radiotracer uptake to neighbours of cells, and when multiplexed with additional radiotracers, fluorescent probes, or in situ tissue analysis, autoradiography can help to characterize the mechanism of radiotracer uptake and assess functional heterogeneity in tissue. In this chapter, the author outlines the basic ex vivo autoradiography protocol and shows how it can be multiplexed using dual radionuclides 18F and 14C. They also highlight where autoradiography can be combined with other technologies to provide synergistic information for interrogating spatial biology.

Beyond Average: α-Particle Distribution and Dose Heterogeneity in Bone Metastatic Prostate Cancer.

α-particle emitters are emerging as a potent modality for disseminated cancer therapy because of their high linear energy transfer and localized absorbed dose profile. Despite great interest and pharmaceutical development, there is scant information on the distribution of these agents at the scale of the α-particle pathlength. We sought to determine the distribution of clinically approved [223Ra]RaCl2 in bone metastatic castration-resistant prostate cancer at this resolution, for the first time to our knowledge, to inform activity distribution and dose at the near-cell scale. Methods: Biopsy specimens and blood were collected from 7 patients 24 h after administration. 223Ra activity in each sample was recorded, and the microstructure of biopsy specimens was analyzed by micro-CT. Quantitative autoradiography and histopathology were segmented and registered with an automated procedure. Activity distributions by tissue compartment and dosimetry calculations based on the MIRD formalism were performed. Results: We revealed the activity distribution differences across and within patient samples at the macro- and microscopic scales. Microdistribution analysis confirmed localized high-activity regions in a background of low-activity tissue. We evaluated heterogeneous α-particle emission distribution concentrated at bone-tissue interfaces and calculated spatially nonuniform absorbed-dose profiles. Conclusion: Primary patient data of radiopharmaceutical therapy distribution at the small scale revealed that 223Ra uptake is nonuniform. Dose estimates present both opportunities and challenges to enhance patient outcomes and are a first step toward personalized treatment approaches and improved understanding of α-particle radiopharmaceutical therapies.

From nuclear track characterization to machine learning based image classification in neutron autoradiography for boron neutron capture therapy.

Knowledge of the 10B microdistribution is of great relevance in BNCT studies. Since 10B concentration assesment through neutron autoradiography depends on the correct quantification of tracks in a nuclear track detector, image acquisition and processing conditions should be controlled and verified, in order to obtain accurate results to be applied in the frame of BNCT. With this aim, an image verification process was proposed, based on parameters extracted from the quantified nuclear tracks. Track characterization was performed by selecting a set of morphological and pixel-intensity uniformity parameters from the quantified objects (area, diameter, roundness, aspect ratio, heterogeneity and clumpiness). Their distributions were studied, leading to the observation of varying behaviours in images generated by different samples and acquisition conditions. The distributions corresponding to samples coming from the BNC reaction showed similar attributes in each analyzed parameter, proving to be robust to the experimental process, but sensitive to light and focus conditions. Considering those observations, a manual feature extraction was performed as a pre-processing step. A Support Vector Machine (SVM) and a fully dense Neural Network (NN) were optimized, trained, and tested. The final performance metrics were similar for both models: 93%-93% for the SVM, vs 94%-95% for the NN in accuracy and precision respectively. Based on the distribution of the predicted class probabilities, the latter had a better capacity to reject inadequate images, so the NN was selected to perform the image verification step prior to quantification. The trained NN was able to correctly classify the images regardless of their track density. The exhaustive characterization of the nuclear tracks provided new knowledge related to the autoradiographic images generation. The inclusion of machine learning in the analysis workflow proves to optimize the boron determination process and paves the way for further applications in the field of boron imaging.

Tumor Control Probability and Small-Scale Monte Carlo Dosimetry: Effects of Heterogenous Intratumoral Activity Distribution in Radiopharmaceutical Therapy.

In radiopharmaceutical therapy, intratumoral uptake of radioactivity usually leads to heterogeneous absorbed dose distribution. The likelihood of treatment success can be estimated with the tumor control probability (TCP), which requires accurate dosimetry, estimating the absorbed dose rate per unit activity to individual tumor cells. Methods: Xenograft cryosections of the prostate cancer cell line LNCaP treated with [177Lu]Lu-PSMA-617 were evaluated with digital autoradiography and stained with hematoxylin and eosin. The digital autoradiography images were used to define the source in a Monte Carlo simulation of the absorbed dose, and the stained sections were used to detect the position of cell nuclei to relate the intratumoral absorbed dose heterogeneity to the cell density. Simulations were performed for 225Ac, 177Lu, and 90Y. TCP was calculated to estimate the mean necessary injected activity for a high TCP. A hypothetical case of activity mainly taken up on the tumor borders was generated and used to simulate the absorbed dose. Results: The absorbed dose per decay to tumor cells was calculated from the staining and simulation results to avoid underestimating the tumor response from low absorbed doses in tumor regions with low cell density. The mean of necessary injected activity to reach a 90% TCP for 225Ac, 177Lu, and 90Y was found to be 18.3 kBq (range, 18-22 kBq), 24.3 MBq (range, 20-29 MBq), and 5.6 MBq (range, 5-6 MBq), respectively. Conclusion: To account for the heterogeneous absorbed dose generated from nonuniform intratumoral activity uptake, dosimetry models can estimate the mean necessary activity to reach a sufficient TCP for treatment response. This approach is necessary to accurately evaluate the efficacy of suggested radiopharmaceuticals for therapy.

Gradients of neurotransmitter receptor expression in the macaque cortex.

Dynamics and functions of neural circuits depend on interactions mediated by receptors. Therefore, a comprehensive map of receptor organization across cortical regions is needed. In this study, we used in vitro receptor autoradiography to measure the density of 14 neurotransmitter receptor types in 109 areas of macaque cortex. We integrated the receptor data with anatomical, genetic and functional connectivity data into a common cortical space. We uncovered a principal gradient of receptor expression per neuron. This aligns with the cortical hierarchy from sensory cortex to higher cognitive areas. A second gradient, driven by serotonin 5-HT1A receptors, peaks in the anterior cingulate, default mode and salience networks. We found a similar pattern of 5-HT1A expression in the human brain. Thus, the macaque may be a promising translational model of serotonergic processing and disorders. The receptor gradients may enable rapid, reliable information processing in sensory cortical areas and slow, flexible integration in higher cognitive areas.

Measurement of an evaporation coefficient in tissue sections as a correction factor for 10B determination.

Boron neutron capture therapy (BNCT) is a cancer treatment option that combines preferential uptake of a boron compound in tumors and irradiation with thermal neutrons. For treatment planning, the boron concentration in different tissues must be considered. Neutron autoradiography using nuclear track detectors (NTD) can be applied to study both the concentration and microdistribution of boron in tissue samples. Histological sections are obtained from frozen tissue by cryosectioning. When the samples reach room temperature, they undergo an evaporation process, which leads to an increase in the boron concentration. To take this effect into account, certain correction factors (evaporation coefficients, CEv) must be applied. With this aim, a protocol was established to register and analyze mass variation of tissue sections, measured with a semimicro scale. Values of ambient temperature, pressure, and humidity were simultaneously recorded. Reproducible results of evaporation curves and CEv values were obtained for different tissue samples, which allowed the systematization of the procedure. This study could contribute to a more precise determination of boron concentration in tissue samples through the neutron autoradiography technique, which is of great relevance to make dosimetric calculations in BNCT.

Synthesis and evaluation of 111 In-labeled tetrapeptide-based compounds as single-photon emission computed tomography imaging probes targeting granzyme B.

Granzyme B is an attractive target as a biomarker for contributing to improve the treatment with immune checkpoint inhibitor (ICI). In this study, we designed novel 111 In-labeled granzyme B-targeting single-photon emission computed tomography (SPECT) imaging probes, [111 In]IDT and [111 In]IDAT. Nonradioactive In-labeled granzyme B-targeting compounds ([nat In]IDT, [nat In]IDAT) showed the affinity for recombinant mouse granzyme B. [111 In]IDT and [111 In]IDAT were obtained with moderate radiochemical yield and high stability in mouse plasma (>95%). In a biodistribution experiment using tumor-bearing mice, [111 In]IDT and [111 In]IDAT showed moderate accumulation in tumor. Ex vivo autoradiography (ARG) indicated that the accumulation of radioactivity in tumor was correlated to expression of granzyme B confirmed by the immunohistochemical staining. These results indicated that [111 In]IDT and [111 In]IDAT showed the basic properties as granzyme B-targeting SPECT probes.

Yttrium-90 Activity Quantification in PET/CT-Guided Biopsy Specimens from Colorectal Hepatic Metastases Immediately after Transarterial Radioembolization Using Micro-CT and Autoradiography.

PURPOSE: To evaluate the yttrium-90 (90Y) activity distribution in biopsy tissue samples of the treated liver to quantify the dose with higher spatial resolution than positron emission tomography (PET) for accurate investigation of correlations with microscopic biological effects and to evaluate the radiation safety of this procedure. MATERIALS AND METHODS: Eighty-six core biopsy specimens were obtained from 18 colorectal liver metastases (CLMs) immediately after 90Y transarterial radioembolization (TARE) with either resin or glass microspheres using real-time 90Y PET/CT guidance in 17 patients. A high-resolution micro-computed tomography (micro-CT) scanner was used to image the microspheres in part of the specimens and allow quantification of 90Y activity directly or by calibrating autoradiography (ARG) images. The mean doses to the specimens were derived from the measured specimens' activity concentrations and from the PET/CT scan at the location of the biopsy needle tip for all cases. Staff exposures were monitored. RESULTS: The mean measured 90Y activity concentration in the CLM specimens at time of infusion was 2.4 ± 4.0 MBq/mL. The biopsies revealed higher activity heterogeneity than PET. Radiation exposure to the interventional radiologists during post-TARE biopsy procedures was minimal. CONCLUSIONS: Counting the microspheres and measuring the activity in biopsy specimens obtained after TARE are safe and feasible and can be used to determine the administered activity and its distribution in the treated and biopsied liver tissue with high spatial resolution. Complementing 90Y PET/CT imaging with this approach promises to yield more accurate direct correlation of histopathological changes and absorbed dose in the examined specimens.

Relationships between neurotransmitter receptor densities and expression levels of their corresponding genes in the human hippocampus.

Neurotransmitter receptors are key molecules in signal transmission, their alterations are associated with brain dysfunction. Relationships between receptors and their corresponding genes are poorly understood, especially in humans. We combined in vitro receptor autoradiography and RNA sequencing to quantify, in the same tissue samples (7 subjects), the densities of 14 receptors and expression levels of their corresponding 43 genes in the Cornu Ammonis (CA) and dentate gyrus (DG) of human hippocampus. Significant differences in receptor densities between both structures were found only for metabotropic receptors, whereas significant differences in RNA expression levels mostly pertained ionotropic receptors. Receptor fingerprints of CA and DG differ in shapes but have similar sizes; the opposite holds true for their "RNA fingerprints", which represent the expression levels of multiple genes in a single area. In addition, the correlation coefficients between receptor densities and corresponding gene expression levels vary widely and the mean correlation strength was weak-to-moderate. Our results suggest that receptor densities are not only controlled by corresponding RNA expression levels, but also by multiple regionally specific post-translational factors.

Non-invasive regional cerebral blood flow quantification in the 123I-IMP autoradiography using artificial neural network.

PURPOSE: Regional cerebral blood flow (rCBF) quantification using 123I-N-isopropyl-p-iodoamphetamine (123I-IMP) requires an invasive, one-time-only arterial blood sampling for measuring the 123I-IMP arterial blood radioactivity concentration (Ca10). The purpose of this study was to estimate Ca10 by machine learning (ML) using artificial neural network (ANN) regression analysis and consequently calculating rCBF and cerebral vascular reactivity (CVR) in the dual-table autoradiography (DTARG) method. MATERIALS AND METHODS: This retrospective study included 294 patients who underwent rCBF measurements through the 123I-IMP DTARG. In the ML, the objective variable was defined by the measured Ca10, whereas the explanatory variables included 28 numeric parameters, such as patient characteristic values, total injection 123I-IMP radiation dose, cross-calibration factor, and the distribution of 123I-IMP count in the first scan. ML was performed with training (n = 235) and testing (n = 59) sets. Ca10 was estimated in testing set by our proposing model. Alternatively, the estimated Ca10 was also calculated via the conventional method. Subsequently, rCBF and CVR were calculated using estimated Ca10. Pearson's correlation coefficient (r-value) for the goodness of fit and the Bland-Altman analysis for assessing the potential agreement and bias were performed between the measured and estimated values. RESULTS: The r-value of Ca10 estimated by our proposed model was higher compared with the conventional method (0.81 and 0.66, respectively). In the Bland-Altman analysis, mean differences of 4.7 (95% limits of agreement (LoA): -18-27) and 4.1 (95% LoA: -35-43) were observed using proposed model and the conventional method, respectively. The r-values of rCBF at rest, rCBF after the acetazolamide challenge, and CVR calculated using the Ca10 estimated by our proposed model were 0.83, 0.80 and 0.95, respectively. CONCLUSION: Our proposed ANN-based model could accurately estimate the Ca10, rCBF, and CVR in DTARG. These results would enable non-invasive rCBF quantification in DTARG.

Assay of GTPγS Binding in Autoradiography.

Autoradiography of radiolabeled GTPγS ([35S]GTPγS) binding is a relevant technique to study the function of G protein-coupled receptors (GPCRs) ex vivo. Here, we describe the protocol for such a method, suitable for investigating CB1 receptor functionality in tissue slices from rodent brains.

Spatial quantification of single cell mRNA and ligand binding of the kappa opioid receptor in the rat hypothalamus.

Detailed quantification of brain tissue provides a deeper understanding of changes in expression and function. We have created a pipeline to study the detailed expression patterns of the kappa opioid receptor in the rat hypothalamus using high resolution fluorescence microscopy and receptor autoradiography. The workflow involved structured serial sampling of rat hypothalamic nuclei, in situ detection of mRNA and receptor expression, and advanced image analysis. Our results demonstrate how maintaining spatial information can lead to increased understanding of RNA and protein expression. In addition, we show the detailed expression patterns of the kappa opioid receptor in the rat hypothalamus.

18F-PSMA Cerenkov Luminescence and Flexible Autoradiography Imaging in a Prostate Cancer Mouse Model and First Results of a Radical Prostatectomy Feasibility Study in Men.

Intraoperative identification of positive resection margins (PRMs) in high-risk prostate cancer (PC) needs improvement. Cerenkov luminescence imaging (CLI) with 68Ga-PSMA-11 is promising, although limited by low residual activity and artificial signals. Here, we aimed to assess the value of CLI and flexible autoradiography (FAR) with 18F-PSMA-1007. Methods: Mice bearing subcutaneous PSMA-avid RM1-PGLS tumors were administered 18F-PSMA-1007, and PET/CT was performed. After the animals had been killed, organs were excised and measured signals in CLI and FAR CLI were correlated with tracer activity concentrations (ACs) obtained from PET/CT. For clinical assessment, 7 high-risk PC patients underwent radical prostatectomy immediately after preoperative 18F-PSMA PET/CT. Contrast-to-noise ratios (CNRs) were calculated for both imaging modalities in intact specimens and after incision above the index lesion. Results: In the heterotopic in vivo mouse model (n = 5), CLI did not detect any lesion. FAR CLI detected a distinct signal in all mice, with a lowest AC of 7.25 kBq/mL (CNR, 5.48). After incision above the index lesion of the prostate specimen, no increased signal was observed at the cancer area in CLI. In contrast, using FAR CLI, a signal was detectable in 6 of 7 patients. The AC in the missed index lesion was 1.85 kBq/mL, resulting in a detection limit of at least 2.06 kBq/mL. Histopathology demonstrated 2 PRMs, neither of which was predicted by CLI or FAR CLI. Conclusion: 18F-PSMA FAR CLI was superior to CLI in tracer-related signal detectability. PC was could be visualized in radical prostatectomy down to 2.06 kBq/mL. However, the detection of PRMs was limited. Direct anatomic correlation of FAR CLI is challenging because of the scintillator overlay.

Small-scale (sub-organ and cellular level) alpha-particle dosimetry methods using an iQID digital autoradiography imaging system.

Targeted radiopharmaceutical therapy with alpha-particle emitters (αRPT) is advantageous in cancer treatment because the short range and high local energy deposition of alpha particles enable precise radiation delivery and efficient tumor cell killing. However, these properties create sub-organ dose deposition effects that are not easily characterized by direct gamma-ray imaging (PET or SPECT). We present a computational procedure to determine the spatial distribution of absorbed dose from alpha-emitting radionuclides in tissues using digital autoradiography activity images from an ionizing-radiation quantum imaging detector (iQID). Data from 211At-radioimmunotherapy studies for allogeneic hematopoietic cell transplantation in a canine model were used to develop these methods. Nine healthy canines were treated with 16.9-30.9 MBq 211At/mg monoclonal antibodies (mAb). Lymph node biopsies from early (2-5 h) and late (19-20 h) time points (16 total) were obtained, with 10-20 consecutive 12-µm cryosections extracted from each and imaged with an iQID device. iQID spatial activity images were registered within a 3D volume for dose-point-kernel convolution, producing dose-rate maps. The accumulated absorbed doses for high- and low-rate regions were 9 ± 4 Gy and 1.2 ± 0.8 Gy from separate dose-rate curves, respectively. We further assess uptake uniformity, co-registration with histological pathology, and requisite slice numbers to improve microscale characterization of absorbed dose inhomogeneities in αRPT.

Correspondence between gene expression and neurotransmitter receptor and transporter density in the human brain.

Neurotransmitter receptors modulate signaling between neurons. Thus, neurotransmitter receptors and transporters play a key role in shaping brain function. Due to the lack of comprehensive neurotransmitter receptor/transporter density datasets, microarray gene expression measuring mRNA transcripts is often used as a proxy for receptor densities. In the present report, we comprehensively test the spatial correlation between gene expression and protein density for a total of 27 neurotransmitter receptors, receptor binding-sites, and transporters across 9 different neurotransmitter systems, using both PET and autoradiography radioligand-based imaging modalities. We find poor spatial correspondences between gene expression and density for all neurotransmitter receptors and transporters except four single-protein metabotropic receptors (5-HT1A, CB1, D2, and MOR). These expression-density associations are related to gene differential stability and can vary between cortical and subcortical structures. Altogether, we recommend using direct measures of receptor and transporter density when relating neurotransmitter systems to brain structure and function.

Autoradiography on deparaffinized tissue sections - A feasibility study with 68Ga-labeled PET-tracers.

Tissue available for retrospective research questions is often already paraffin-embedded for better preservation. However, in vitro autoradiography (AURA) is normally performed on cryopreserved tissue sections. We hypothesized a) that it would also be feasible with deparaffinized tissue sections, enabling the use of human paraffin-embedded tissue for in vitro AURA and b) that the results would be comparable to those obtained with corresponding cryosections. For that purpose, the clinically relevant oncological targets CXCR4, SSTR and PSMA were evaluated. In vitro AURA on deparaffinized tissue sections was feasible, but only with the two receptor ligands [68Ga]Ga-PentixaFor and [68Ga]Ga-DOTANOC. [68Ga]Ga-PSMA-11 did not show any binding on deparaffinized tissue sections, suggesting that native tissue is required for an interaction between this inhibitor and the enzyme.

Quantifying 226Ra activity in a complex assemblage of 226Ra-bearing minerals using alpha autoradiography and SEM/EDS.

226Ra is an ultra-trace element with important environmental implications for many industries (including water treatment and oil and mineral extraction). Its extremely low concentrations in natural environments do not allow for direct observation and measurement of the 226Ra-bearing minerals governing 226Ra mobility. To better understand the retention processes for 226Ra in rocks and soil, a synthesized assemblage of 226Ra-doped minerals was made, combining montmorillonite, ferrihydrite and barite. A new methodology was developed using alpha activity maps acquired using alpha autoradiography, and elemental maps by using SEM/EDS. These maps were processed using a global approach, considering the entirety of the signal. The comparison of the alpha activity map and the elemental map enabled a correlation to be established between the 226Ra activity and the chemical composition and identification of the main 226Ra-bearing mineral of the assemblage, from which we were able to estimate the contribution of each mineral to the total activity of the assemblage, and to quantify the 226Ra-activity for each mineral. This methodology makes it possible to link mineralogy and occurrence of 226Ra at the scale of the mineral (tens of µm). It can be applied to natural samples, including fine-grained samples with a complex mineralogy.

Orexin 2 receptor (OX2R) protein distribution measured by autoradiography using radiolabeled OX2R-selective antagonist EMPA in rodent brain and peripheral tissues.

Orexin, a neuropeptide, performs various physiological functions, including the regulation of emotion, feeding, metabolism, respiration, and sleep/wakefulness, by activating the orexin 1 receptor and orexin 2 receptor (OX2R). Owing to the pivotal role of OX2R in wakefulness and other biological functions, OX2R agonists are being developed. A detailed understanding of OX2R protein distribution is essential for determining the mechanisms of action of OX2R agonists; however, this has been hindered by the lack of selective antibodies. In this study, we first confirmed the OX2R-selective binding of [3H]-EMPA in in vitro autoradiography studies, using brain slices from OX2R knockout mice and their wild-type littermates. Subsequently, OX2R protein distribution in rats was comprehensively assessed in 51 brain regions and 10 peripheral tissues using in vitro autoradiography with [3H]-EMPA. The widespread distribution of OX2R protein, including that in previously unrecognized regions of the retrosplenial cortex, was identified. In contrast, OX2R protein expression was negligible/very low in peripheral tissues, suggesting that orexin exerts OX2R-dependent physiological functions primarily through activation of the central nervous system. These findings will be useful for understanding the wide range of biological functions of OX2R and the application of OX2R agonists in various disorders.